AD-A178  761  SIGNATURES  FROM  IN-SITU  LIGHT  SCATTERING  OF  1/1 

'  NON-HOHOGENEOUS  NON-SPHERICAL  (U)  GEORGIA  INST  OF  TECH 

ATLANTA  SCHOOL  OF  GEOPHYSICAL  SCIENCES  A  COLETTI 
UNCLASSIFIED  0b  FEB  87  ARO-21016  1-GS  DAAG29-84-K-0115  F/G  20/6  NL 


fl-dc  J/O/U-QS 


FINAL  REPORT 


m  Fill  Cob- 


signatures  FROM  IN-SITU  LIGHT  SCATTERING  OF 
NON-HOMOGENEOUS  NON-SPHERICAL  PARTICLES 


By 

Alessandro  Coletti 


February  6,  1 987 
Prepared  for 

U.S.  ARMY  RESEARCH  OFFICE 

Under 

Contract  Number  DAAG29-84-K>01 1 5 

GEORGIA  INSTITUTE  OF  TECHNOLOGY 

A  UNIT  OF  THE  UNIVERSITY  SYSTEM  OF  GEORGIA 
SCHOOL  OF  GEOPHYSICAL  SCIENCES 
ATLANTA,  GEORGIA  30332 


$7  4  l  2$3 

APPROVED  FOR  PUBLIC  RELEASE; 

DISTRIBUTION  UNLIMITED 


OTIC 

iELECTEI 


SIGNATURES  FROM  IN-SITU  LIGHT  SCATTERING  OF 
NON-HOMOGENEOUS  NON-SPHERICAL  PARTICLES 


FINAL  REPORT 


Alessandro  Coletti 


February  6,  1 987 


U.S.  ARMY  RESEARCH  OFFICE 


Contract  Number  DAAG29-84-K-01 1 5  \ 


SCHOOL  OF  GEOPHYSICAL  SCIENCES 
GEORGIA  INSTITUTE  OF  TECHNOLOGY 
ATLANTA,  GEORGIA  30332 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


SCCUMITV  CLASSIFICATION  OF  THIS  FACE  (Wbon  O—  Bntmnd) 


REPORT  DOCUMENTATION  PAGE 


2.  OOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


}  RECIPIENT'S  CATALOG  NUMBER 


A.  TITLE  (*>d  Svbllllo) 


Signatures  from  in-si tu  light  scattering  of 
non-homogeneous  non  spherical  particles. 


PFRFOMMIM6  QRG  REPORT  NUMBER 


7.  AU  T  HORf  9) 


Alessandro  S.  Coletti 


•.PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

School  of  Geophysical  Sciences 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 


M.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

U.  S.  Army  Research  Office 
Post  Office  Box  12211 

Research  Triangle  Park,  NC  27709 _ 


.  MONITORING  AGENCY  NAME  •  ADORES VII  dllloronl  (torn  Controlling  Oil! CO) 


a  contract  or  grant  numberoi 

DAAG  29-84-k-01 15 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  a  WORK  UNIT  NUMBERS 


I 


IS.  NUMBER  OF  PACES 


IS.  SECURITY  CLASS,  (ol  IW«  report) 


Unclassified 


is*,  oecl assification/downgraoing 
schedule 


IS.  DISTRIBUTION  statement  (ol  (Ala  Ropori; 


Approved  for  public  release;  distribution  unlimited. 


17  DISTRIBUTION  STATEMENT  (ol  iho  obotroct  nlwM  In  Block  to,  II  dllloronl  trmmi  Rapa rt) 


It.  SUPPLEMENTARY  notes 

The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
,  position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 

tf,_  KEY  vobOS  fCwiUnu*  on  footrot  old#  it  nocoo«»f  lamnttty  oy  oiocm 

—»  Light  scattering,*  aerosol  particles',  laser  measurements,' numerical  inversions. 

r,  i  e  n/,*ri/'\i;trP  ^ 

\ 


jj|  iwmirt  pTwir^ii  mm  1 1 - rrr~  M  '  timMlty 

The  research  performed  under  this  program  has  focused  on  studying  the  angular 
light  scattering  of  monodisperse  quasi  spherical  non  homogeneous  particles.  For 
that  purpose  we-had  to  develop  specific  laboratory  techniques ^suitable  for  the 
generation  of  these  particles  and  for  properly  collecting  the  data.  For 
comparing  the  measurements  with  the  light  scattering  theories  for  spheres  and 
spheroids  we  have  formulated  a  new  kind  of  algorithm  that  is  innovative  in  many 
aspects  and  has  useful  application.  The  method  might  be  advantageous  for, 


00  /jSTt »  1473  EDITION  OF 


•  MOV  ••  IS  OBSOLETE 


UNCLASSIFIED _ _  /  1 

SECURITY  Cl  AMI  tic  at  ion  of  This  PAGE  fWNon  P«l«  Bnlorod) 


UNCLASSIFIED 


\ 


(20.  ABSTRACT  cont.) 


inversions  of  physical  parameters  from  experimental  data,  and  for 
estimating  '  their  statistical  significance.  ■'  li. 

*■*-  (a  /  '  ^ 


f* 


security  classification  Of  tmis  AA5(m«i  o»*»  *"*•»•«) 


TABLE  OF  CONTENTS 


Program  Overview  .  1 

1.  Introduction  .  3 

2.  Applications  .  8 

2.1  The  Direct  Problem . 8 

2.2  The  Inverse  Problem . 10 

3.  Conclusions . . . 15 

4.  List  of  Technical  Reports  Published . 16 

5.  List  of  Participating  Personnel . 16 

References . 17 


Appendix  A 

Aerosol  Composition  (8  pages) 


Appendix  B 

Fisheye  Polar  Nephelometer  (11  pages) 


Appendix  C 

Measured  Phase  Functions  (24  pages) 


it/lJ 


Accession  For 

1 

n 


By _ 

Dirt- 


i  -  t 


itti  CJ 

rTIC  TV-; 
Un.f.m.v 
Ju.;t  i . 


PROGRAM  OVERVIEW 


The  major  objective  of  this  research  is  to  further  develop 
numerical  methods  for  the  analytical  treatment  of  angular  light 
scattering  from  randomly  oriented  non-spherical  non-homogeneous 
particles.  This  requirement  has  been  achieved  with  the 
development  of  the  numerical  method  described  in  the  following 
paragraphs  which  is  the  result  of  our  study  performed  on 
laboratory  generated  light  scattering  data.  We  feel  that  the 
approach  we  are  suggesting  has  many  innovative  aspects,  and  can 
be  used  in  many  practical  problems. 

During  the  work  we  have  devised  laboratory  techniques 
(described  in  Appendix  A)  for  generating  the  data,  and  we  have 
performed  a  feasibility  study  on  a  new  kind  of  light  scattering 
instruments  (described  in  the  Appendix  B)  that  performs  angular 
scans  without  using  moving  optical  parts.  Homogeneous  spherical 
particles  made  with  mixed  solutions  of  ammonium  sulfate- 
methylene  blue  and  ammonium  sulfate-  nigrosine  dye  have  been  used 
as  a  reference.  Inhomogeneous  particles  have  been  made  by  mixing 
ionic  salts  with  gold  and  silver  colloids  at  various  stages  of 
their  nucleation.  Metallic  colloids  of  silver  and  gold  have  been 
chosen  because  they  nucleate  in  relatively  narrow  size 
distributions  and  relatively  small  uncertainties  exist  in  the 
values  of  the  refractive  index.  Some  typical  results  of  angular 
light  scattering  measurements  obtained  with  these  chemical 
solutions  are  reported  in  Appendix  C  of  this  note. 
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We  have  confirmed  that  the  angular  light  scattering 
measurements  taken  every  5°  in  the  angular  interval  between  10 
and  170°,  could  not  be  interpreted  with  the  Lorenz-Mie  theory  for 
homogeneous  spherical  particles.  However,  the  application  of 
light  scattering  theory  for  spheroids  (with  the  EBCM  program 
kindly  provided  to  us  by  Dr.  Peter  Barber)  turned  out  to  be 
cumbersome  and  of  doubtful  use.  In  fact  in  the  few  cases  where 
the  particle  sizes  were  small  enough  to  permit  practical 
applications  of  the  EBCM  theory,  limitations  existed  on  the 
statistical  significance  of  the  particle  sizes  and  axial  ratios 
as  obtained  with  the  electron  scanning  microscope,  and  on  the 
values  of  the  effective  refractive  index  applicable  to  the 
mixtures . 
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Consequently,  a  different  approach  has  been  attempted  where 
the  optical  measurements  are  analyzed  for  their  information 
content  in  terms  of  particle  composition  and  shape.  Basically 
the  method  consists  of  calculating  the  Legendre  polynomial  series 
fitting  the  experimental  data.  The  set  of  these  coefficients  can 
be  compared  with  theoretical  sets  calculated  for  different 
properties  of  the  particles.  The  advantage  of  the  method, 
against  standard  inversion  methods,  arises  mostly  from  the  more 
compact  representation  of  the  phase  functions  allowed  by  the  use 
of  the  Legendre  polynomials.  In  what  follows  we  report  a  few 
examples  of  how  the  proposed  method  may  be  applied  in  practice. 


►  '  >  V \ 


1.  INTRODUCTION 


Angular  light  scattering  measurements  of  natural  and 
laboratory  generated  aerosols  contain  information  on  the  size, 
chemical  composition  and  shape  of  the  particles.  Nevertheless, 
up  to  the  present  date,  the  numerical  inversion  techniques 
applied  to  experimental  data  have  been  based  on  the  Lorenz-Mie 
theory  for  spherical  homogeneous  particles  (Ferrara  et  al.,  1970; 
Grams  et  al.,  1974;  Quenzel,  1970;  Takamura  and  Tanaka,  1985). 
To  the  best  of  our  knowledge,  the  only  numerical  study  on  the 
possibility  of  inferring  particle  shapes  from  simulated  light 
scattering  data  has  been  made  by  Heintzenberg  and  Welch  (1982). 
These  authors  show  that  even  though  the  light  scattering  contains 
information  on  the  particle  shape,  still  its  retrieval  is 
practically  impossible  if  the  assumption  is  made  that  no  previous 
knowledge  of  particle  shape  is  available. 

The  differences  in  the  light  scattering  by  spheres  and 
spheroids  have  been  described  by  Yamamoto  and  Tanaka  (1969)  and 
Barber  and  Yeh  (1975).  They  have  shown  that  basic  quantitative 
and  qualitative  differences  exist  even  though  the  position  of 
maxima  and  minima  in  both  scattering  coefficient  and  asymmetry 
parameter  (proportional  through  a  numerical  constant  to  the  first 
coefficient  of  the  Legendre  polynomial  series)  for  spheroidal  and 
spherical  particles  have  some  qualitative  similarities. 
Important  contributions  to  the  numerical  capability  for  computing 
the  light  scattering  properties  of  non-spher ical  particles  have 


been  achieved  by  Barber  and  Liou  (see  e.g.  Liou  et  al.  1983). 
They  also  focused  their  attention  on  some  of  the  particle 


characteristics  that  can  facilitate  their  recognition  in 
experiments . 

Basically,  the  data  analysis  described  here  provides 
criteria  for  the  recognition  of  non-spherical  particles  from 
light  scattering  data.  The  method  follows  from  calculations  made 
by  Chu  and  Churchill  (1955)  on  the  Legendre  polynomial  expansion 
of  the  Lorenz  Mie  theory  and  discussions  by  Van  de  Hulst  (1980). 
Interest  in  the  study  of  the  Legendre  polynomial  coefficients 
originates  both  from  their  advantageous  mathematical  properties 
and  from  their  physical  meaning  in  the  spherical  harmonics 
expansion  of  the  scattered  field. 

Fig.  1  shows  some  examples  of  Legendre  polynomial 
coefficients  computed  for  spherical  particles  of  various 
refractive  indices  as  a  function  of  one  half  of  their  phase  shift 
p  (p=4wr/*  [m-1]).  The  various  curves  represent  the  coefficients 
from  the  1st  to  the  19th  order  polynomial.  Fig.  1  a)  is  for  a 
refractive  index,  m=1.3,  b)  for  m=1.4,  c)  for  m=1.5,  and  d)  for  a 
light  absorbing  particle  with  m=1.5  -  i0.5.  The  compactness  of 
this  kind  of  representation  can  be  easily  appreciated  when  it  is 
realized  that  each  of  these  figures  describes  all  the  major 
features  of  the  phase  functions  in  the  interval  under 
consideration.  The  resolution  used  in  the  figures  shows  some  of 
the  characteristics  of  the  Legendre  polynomial  coefficients  that 
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Van  de  Hulst  could  not  determine  from  Chu  and  Churchill 


calculations.  Namely,  the  maxima  and  minima  of  the  coefficients 
appear  to  occur  at  similar  values  of  phase  shift  until  the  order 
of  the  coefficient  is  larger  than  x-1/2  (x=2nr/x);  for  higher 
values  of  the  order,  their  position  moves  with  an  almost  power 
law  toward  larger  values  of  p  as  the  order  of  the  coefficient 
increases . 

Numerically,  the  Legendre  polynomial  coefficients  of  both 
spherical  and  spheroidal  particles  can  be  easily  obtained  from 
the  existing  numerical  codes  since  the  Legendre  polynomials  are 


readily  available  as  part  of  these  of  programs.  They  can  instead 
be  estimated  experimentally  from  angular  light  scattering 
measurements,  but  only  under  specific  experimental  conditions. 


For  a  given  phase  function,  4>(p),  the  n-th  order  coefficient,  un. 


is  given  by: 
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(2n+l)  to 


4> 

.  -1 


(1 


where  PR(p)  is  the  n-th  order  Legendre  polynomial  and  p  =  cos  0. 


If  we  have  N  light  scattering  measurements,  (i  =  l,..,N), 


performed  at  angles  0^ ,  we  can  compute  N  parameters  gR : 
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Z.  V.w.P  (P. ) 
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measures  light  scattering  phase  functions  within  a  given  angular 
interval.  For  larger  particle  sizes  where  the  contributions  due 
to  polynomials  of  higher  order  become  sensible,  the  parameters  g^ 
measured  in  a  limited  interval  of  angles  can  not  be  any  longer 

Therefore  we  can  compute  only  the 
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where  and  are  the  cosines  of  the  minimum  and  maximum 
measurable  angles  and  R  the  remainder.  It  will  be  shown  in  the 
next  section  that  even  in  this  case  there  is  a  usable  information 
content  of  the  g^  coefficients  of  the  various  orders  so  that, 
under  suitable  experimental  conditions,  it  is  still  possible  to 
overcome  some  of  the  drawbacks  described  by  Nugnai  and  Wiscombe 
and  by  Heintzenberg  and  Welch. 

In  the  meantime,  we  can  observe  that  independently  of 
whether  eq.  2  or  eq.  5  are  satisfied,  the  optimized  Legendre 
polynomial  series  of  experimental  data  provides  the  best 
interpolation  method  of  the  measurements  independent  of  particle 
shapes  in  the  least  square  sense.  This  interpolation  is 
preferable  to  the  one  we  have  previously  proposed  (Coletti,  1984) 
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where  interpolation  of  the  forward  light  scattering  of 
experimental  phase  function  of  quasi  spherical  particles  was 
performed  using  diffraction  theory,  and  where  the  backward 
scattering,  important  for  remote  sensing  measurements,  was 
treated  empirically. 

2.  APPLICATIONS 

2.1  The  direct  problem 

To  illustrate  the  practical  advantages  presented  by  the 
representations  in  terms  of  gQ  coefficients,  let  us  consider  the 
following  problem:  we  want  to  perform  angular  light  scattering 
measurements  to  distinguish  between  monodisperse  spherical 
particles  with  real  refractive  index  1.5  from  particles  with 
complex  refractive  index  1.5  -iO.Ol  in  the  size  parameter  range 


between  3  and  20.  We  want  to  estimate  the  best  angles  to  use  for 
measurements  where  no  dependable  measurements  can  be  taken  at 
angles  smaller  than  10°  and  greater  than  170°. 

For  the  purpose  we  suppose  that  the  relative  experimental 
error,  €,  has  a  statistical  distribution  with  variance 
proportional  to  the  true  values  of  the  V’s,  V*.  Then  we  can  see 
that  (dropping  the  weights  in  the  summation  for  simplicity): 


N  N  N 

E.  v.p  .=  E.  V  .P  .+  eE.V  .  P 
xi  x  n , x  ji  l  n,i  j  l  n,i 
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That  is,  under  our  assumption,  the  relative  errors  in  the  gn  are 
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the  sane  for  any  n.  Therefore  in  Fig.  2,  where  the  relative 

variations  of  the  gn  coefficients  in  the  two  cases  of  the  light 

absorbing  and  non  light  absorbing  particles  are  shown  (for  the 

refractive  indices  1.5  and  1.5  -iO.Ol),  the  two  kinds  of 

particles  will  be  indistinguishable  whenever  the  relative 

differences  of  all  the  considered  coefficients  will  be 

simultaneously  smaller  than  the  errors.  Fig.  2  has  been  drawn  by 

computing  the  g's  for  the  two  kinds  of  particles  in  the  angular 

interval  between  10  and  170°.  The  isopleths  have  only  an 

illustrative  purpose  since  the  functions  gn(x)  do  not  describe  a 

continuous  surface.  Following  the  example  under  discussion,  we 

can  estimate  that  performing  measurements  with  a  precision  of  10* 

we  need  to  estimate  g  coefficients  up  to  the  5th  order  to 

n 

discriminate  light  absorbing  and  non  light  absorbing  particles  in 
the  size  parameter  interval  between  6  and  17.  In  fact,  for  any 
of  those  particles,  Fig  .2  shows  that  there  is  at  least  one  of 
the  first  five  coefficients  that  differs  in  the  two  cases  for 
more  than  the  expected  experimental  error. 

Numerical  procedures,  analogous  to  the  one  just  described, 
can  be  used  to  obtain  the  sensitivity  of  instruments  to  detect 
significant  variations  in  the  particle  properties.  With  the 
traditional  method  quantities  such  as: 
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would  have  been  evaluated  to  estimate  if  the  differences  existing 
between  two  phase  functions,  and  ,  are  sufficiently  large  to 
be  detected.  No  information  would  have  been  obtained  on  the 
differences  existing  between  the  spherical  harmonic  structures  of 
the  phase  functions,  therefore  no  indication  would  become 
available  with  this  latter  method  on  how  to  optimize  the  number 
and  the  position  of  the  angles  for  the  measurements. 

If  needed,  evaluations  similar  to  those  made  in  Fig.  2  can 
be  made  also  in  a  variety  of  other  cases  where  it  is  desired  to 
estimate  the  capability  to  resolve  between  particles  of  different 
shapes,  composition,  aggregation  state,  or  polydispersion. 

2.2  The  Inverse  Problem 

We  now  consider  the  simple  exercise  where  we  want  to 
determine  with  the  Lorenz-Mie  theory  for  spherical  homogeneous 
particles  the  effective  refractive  index  of  a  monodisperse 
aerosol  of  sodium  chloride  particles.  The  numerical  implications 
of  the  method  for  general  inversions  will  be  discussed  later. 

The  Sodium  Chloride  particles  have  a  quas i-spher ical  shape 
and  hollow  internal  structure  as  discussed  by  Coletti  (1982),  the 
light  scattering  phase  function  for  a  particle  of  size  parameter 
x=  5  is  shown  in  Fig.  3.  The  squares  in  the  figure  represent  the 
results  of  the  measurements.  The  continuous  line  in  the  figure 
represent  the  results  of  a  best  fit  obtained  with  the  Lorenz-Mie 
theory  for  a  size  parameter  5  and  a  refractive  index  1.54.  We 
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can  see  in  the  figure  that  the  fit  is  not  totally  satisfactory, 
and  we  want  to  investigate  whether  there  exists  an  effective, 
purely  real,  value  of  refractive  index  that  can  provide  a  better 
fit . 

We  compute  the  optimized  Legendre  polynomial  series  that 
best  fits  the  experimental  values,  and  we  compare  the 
coefficients  with  those  of  the  best  fitting  sphere  in  the  two 
cases  of  refractive  index  1.54  (refractive  index  of  salt)  and 
1.3.  This  latter  value  of  the  refractive  index  has  been  obtained 
applying  the  simple  volume  mixing  of  the  refractive  index  of 
vacuum  and  salt  in  the  ratio  discussed  in  our  previous  work.  Fig 
4  shows  the  coefficient  of  the  optimized  series  as  computed  from 
the  measurements  (squares)  and  the  line  boundaries  determined  by 
the  two  theoretical  spherical  particles  at  the  two  indices  of 
refraction  (upper  continuous  line,  m=1.54;  lower  continuous  line, 
m=1.3).  The  vertical  line  in  the  figure  has  been  drawn  for  the 
value  of  n  for  which  grazing  rays  are  given  by  the  localization 
principle.  We  can  observe  that  the  experimental  coefficients 
have  values  that  shifts  between  the  two  boundaries.  Therefore, 
considering  Fig.  4  and  the  results  discussed  in  Fig  1  a,b,c,  we 
can  conclude  not  only  that  neither  value  of  the  refractive  index 
satisfies  the  requirements  of  a  good  fit,  but  also  that  a  good 
fit  does  not  exist  for  any  other  value  of  refractive  index  in  the 
considered  interval.  Therefore,  the  Lorenz-Mie  theory  for 


I 


spherical  homogeneous  particles  with  x  =  5  and  m=1.54  provides  a 


fit  that  is  dependable  only  up  to  the  1st  order  of  the  polynomial 
expansion  in  the  Legendre  series,  and  the  particles  with  m=1.3 
fit  the  coefficients  with  order  between  7  and  10.  Nevertheless, 
these  features  are  specific  for  sodium  chloride  particles  of  this 
size,  and  it  can  be  seen  in  Fig.  2,  Appendix  C,  that  different 
behavior  is  exhibited  for  other  particle  sizes. 

In  more  a  general  numerical  inversion  problem,  a  minimum 
chi-square  is  generally  searched  assuming  that  each  measured 
value  can  be  represented  as  a  linear  combination  of  a  number 
concentration  of  particles  (having  given  size,  refractive  index, 
shape)  and  their  theoretical  values  of  light  scattering  at  the 
various  angles,  that  is,  using  matrix  notations: 


tf=IIAIl  c 


(8 


where  c  is  the  array  of  the  particle  concentration,  and  the 


elements  a.  .  of  the  matrix  A  are  the  theoretical  values  of  the 
ij 

light  scattering  of  the  particle  of  size  x.  at  the  angle  0.: 

J  ^ 


a  .  .=  ♦(x  . ;  0.  ) 
ij  J  i 


(9 


(10 


The  measurement  array  has,  by  definition,  all  the  Legendre 
polynomial  coefficients  of  order  higher  than  the  total  number  of 


m 


8 
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mm 


the  measurements,  N,  equal  to  zero.  Then,  it  is  easily  seen  that 
this  portion  of  the  series 


£  w  P 


,;b  m  m 


coming  from  the  theoretical  values  of  the  4>’s,  and  normally 


included  in  the  calculations  made  so  far  (see  e.g.  Grams  et  al. 


1974),  can  lead  to  erroneous  results  under  unfavorable 


conditions.  In  fact,  in  the  form  of  eq.  4  this  quantity  escapes 


any  theoretical  analysis  of  aliasing  and  of  the  estimates  of  the 


remainder  of  the  series.  In  contrast,  analysis  with  the  Legendre 


polynomial  along  with  a  suitable  quadrature  formula  (see  e.g. 


Scarborough,  1966)  permits  quantitative  estimates  of  these  latter 


quantities . 


We  can  briefly  discuss  the  direct  use  of  the  coefficients  g 


in  the  case  of  a  general  linear  inversion.  Substituting  eq.  4  in 


eq.  8  we  can  write: 


M©;  ) 


l j  Qn  n  j  n  i 


and  accordingly  introduce  the  matrices  J1  and  ft: 


II A 11  =  II Jill  llftll 


where  the  generic  element  of  the  matrix  ft  is 


i 


wav*1 


and  the  generic  element  of  the  matrix  IT  is: 


T  .  =  P  (0.  ) 

ni  n  1 ' 


(14 


T 

Multiplying  both  sides  of  eq.  8  by  the  matrix  IT  and  by  the 
T  -1 

matrix  HIT  fill  we  obtain: 


iinTnn  1mittii  V  -  unit  c  (15 

But  the  left  hand  side  of  the  equation  is  basically  the  same  of 
eq.  2,  then 

g=  linil  c  (16 

Eq.  16  is  similar  to  eq.  8  and  can  be  solved  with  the  usual 

methods  (S.  Twomey,  1977;  G.  Backus  and  F.  Gilbert,  1970).  The 

use  of  this  equation  instead  of  eq.  8  has  the  further  advantage 

that  the  values  of  the  elements  u.  .  can  be  expected  to  cover  a 

range  of  numerical  values  narrower  than  those  assumed  by  the 

corresponding  elements  A.  ..  It  can  be  anticipated  that  this 

*  J 

should  give  some  numerical  advantage  during  the  calculation  of 

the  inverse  matrix.  In  fact,  the  elements  a.  .  that  represent  the 

*  J 

values  of  the  phase  function  of  a  particle  x.  at  the  angle  6., 

^  J 

often  span  three  or  four  decades,  and  this  makes  cumbersome  the 
already  numerically  delicate  inversion  of  the  matrix  A. 
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The  parameters  obtained  by  the  numerical  inversions,  can  be 
compared  with  the  original  experimental  data  in  a  fashion  totally 
similar  to  that  described  in  the  case  of  monodisperse 
non-homogeneous  particles.  Therefore  similar  insights  into  the 
nature  of  particle  polydispersion  should  be  possible. 

3.  CONCLUSIONS 

The  purpose  of  the  examples  discussed  in  this  note  is  to 
illustrate  some  applications  where  the  series  expansion  of  phase 
functions  in  spherical  harmonics  can  be  applied.  Basically,  we 
have  shown  how  it  is  possible  to  analyze  and  describe  differences 
existing  between  phase  functions,  or  classes  of  phase  functions, 
more  efficiently  then  by  comparing  them  point  by  point.  We  have 
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shown  how  this  can  be  useful  during  analysis  of  light  scattering 
data  of  particle  mixtures  and  to  verify  up  to  what  extent  the 
hypothesis  in  use  properly  approximate  the  physical  or  chemical 
state  of  the  particles. 

We  feel  that  the  use  of  spherical  harmonics  for  studying  the 
light  scattering  properties  of  particles  has  not  been 
sufficiently  explored.  In  fact  the  methods  described  here  are 
just  very  simple  consequences  of  the  fact  that  the  spherical 
harmonics  coefficients  permits  a  more  compact  representation  of 
particle  phase  functions  numerically  and  graphically. 

In  this  paper  objective  criteria  for  optimizing  experimental 
procedures  of  angular  light  scattering  have  been  provided  for  the 
first  time. 
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Fig.  1  -  Legendre  polynoaial  coefficients  for  various  spherical 
particles  as  function  of  one  half  of  the  phase  shift  p  (=2x[m- 
1]):  a)  refractive  index  1.3;  b)  refractive  index  1.4;  c) 
refractive  index  1.5;  c)  refractive  index  1.5  -i  0.5.  Note  in 
the  figure  how  the  minima  occurring  for  the  phase  shifts  near  7 
suddenly  start  shifting  toward  higher  values  when  the  order  of 
the  coefficients  becomes  higher  then  7. 


SIZE  PARAMETER  ,  X 


Fig.  2  -  Line  contours  in  a  table  reporting  relative  variations 
of  the  Legendre  polynomial  coefficients  of  various  orders 
computed  for  different  size  parameters  of  light  absorbing  and  non 
light  absorbing  spherical  particles  (m=1.5  -iO.l  and  1.5). 
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Fig.  4  -  Comparison  of  the  Legendre  polynomial  coefficients  as 
computed  from  the  measurements  in  Fig.  3  and  the  corresponding 
values  obtained  for  spheres  of  size  parameter  5  and  refractive 
index  1.3  (lower  curve)and  1.54  (upper  curve).  The  theoretical 
values  of  the  u  had  been  artificially  connected  with  a 
continuous  line  to  improve  the  clarity  of  the  figure.  The 
vertical  line  for  n=4  indicates  the  order  of  the  grazing  rays  as 
given  by  the  localization  principle. 
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APPENDIX  A 


AEROSOL  COMPOSITION 

To  measure  the  light  scattering  properties  of  non-spherical 
particles  and,  separately,  of  non-homogeneous  particles  we  used 
the  Berglund  Liu  monodisperse  aerosol  generator  that  delivers 
particles  with  an  almost  monodisperse  mass  distributions. 

Non-spherical  particles  have  been  used  to  test  the  light 
scattering  theories  for  spheroidal  particles  (mostly  the  EBCM 
extended  boundary  condition  method  -  for  spheroids)  on 

irregularly  shaped  randomly  oriented  particles.  To  generate 
particles  with  the  desired  characteristics  we  had  to  develop  new 
methods.  In  fact  we  could  not  use  the  usual  methods  of  drying 
droplets  containing  ions  with  strong  polar  bonds  in  solution, 
because  these  particles  often  present  sharp  edges  (e.g.  Leong, 
1985)  and  are  unsuitable  for  the  problem  we  had  to  investigate. 
During  the  measurements  we  have  used  mixtures  of  ionic  salts, 
like  ammonium  sulfate  or  sodium  chloride,  and  organic  dyes,  such 
as  methylene  blue  and  nigrosine  dye.  The  results  we  have 
obtained  using  these  mixtures  indicated  that  these  solutions  when 
evaporated  at  the  appropriate  speed,  yield  smooth  non-spherical 
aerosol  particles.  In  our  case,  the  speed  of  evaporation  was 
controlled  by  varying  the  percentile  of  alcohol  present  in  the 
solvent.  The  microphotograph  reported  in  Fig.  1  shows  as  example 
particles  of  sodium  chloride  and  nigrosine  dye  (mixed  in  equal 
parts  in  a  solvent  50S.-50S:  water  and  methanol).  Fig  2  shows 


anologuous  results  obtained  with  nigrosine  dye  and  aiaoniuB 
sulfate  (mixed  in  a  solvent  75*  methanol-25*  water). 

Quas i-spher ical  particles  containing  homogeneous 
distribution  of  light  absorbing  particles  (with  size  much  smaller 
than  the  wavelength)  have  been  chosen  to  study  the  applicability 
of  the  mixing  rules  for  the  refractive  index.  We  have  tested 
mixtures  of  metallic  colloids  and  ionic  salts  and  we  have 
excluded  particles  containing  carbon  or  soot  particles.  In  fact, 
properly  prepared  Silver  and  Gold  colloids  are  formed  by  narrow 
size  distributions  of  spherical  particles.  The  refractive  index 
of  these  metals  has  been  carefully  studied  and  the  optical 
properties  of  the  colloids  have  been  extensively  described  and 
discussed  in  the  open  literature. 

Gold  colloids  have  been  obtained  using  procedures  described 
by  Turkevich  et  al  (1951).  A  very  narrow  size  distribution  of 
small  particles  is  obtained  from  the  reaction  of  sodium  citrate 
and  chlorauric  acid;  this  colloid  is  subsequently  "developed" 
adding  hydroxylamine  hydrochloride  that  further  coagulates 
metallic  gold  on  the  existing  nuclei.  This  two  step  process 
produces  a  well  defined  narrow  size  distribution  of  spherical 
gold  particles  of  the  desired  size  but  has  a  low  yield  since  it 
is  leaving  a  large  fraction  of  gold  in  ionic  form.  The  reaction 
between  chlorauric  acid  and  citric  acid  has  been  used  for  his 
high  yield  in  gold  particles  even  though  the  particle  size 
distribution  and  shape  were  not  well  defined  as  in  the  previous 
case . 

With  the  reaction  with  the  sodium  citrate,  we  have  obtained 
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a  distribution  of  spheres  centered  at  about  15nm. In  the  reaction 
with  the  citric  acid  a  broad  distribution  of  spheres  and 
triangular  plates  in  the  size  range  within  10-50  nm  has  been 


obtained . 

Silver  colloids  have  been  produced  with  the  reaction  that 
has  been  studied  along  others  by  Von  Raben  et  al.  (1984).  This 
reaction  occurs  at  0°C  between  silver  nitrate  and  sodium 
borohydrade;  the  nucleation  process  has  an  high  yield  in 
spherical  particles  with  narrow  size  distribution  around  25  nm 
diameter  with  a  20%  deviation.  A  good  quality  outcome  from  this 
reaction  critically  depends  on  the  temperature  applied  during  the 
process  and  few  attempts  have  been  necessary  to  produce  a  clear 
red-brownish  sol  of  good  quality  stable  with  time. 

These  metallic  colloids  have  been  used  with  the  vibrating 
orifice  generator  replacing  the  filter  size  suggested  by  the 
manufacturer  (0.05  pm  pores)  with  1  pm  pores  filters.  These 
filters  have  been  found  to  suitably  retain  the  coarser  particles 
while  allowing  the  majority  of  the  particles  with  the  specified 
diameter  to  pass  through.  In  fact  no  significant  change  in  the 
optical  density  of  sample  solutions  could  be  detected  before  and 
after  filtration. 

Unfortunately,  the  analysis  at  the  electron  scanning 
microscope  of  the  aerosol  sampler  collected  on  filters  during  the 
measurements,  gave  negative  results  in  the  case  of  colloidal 
solutions.  In  fact,  these  particles  totally  lost  their  original 
shape  during  the  impact  with  the  filter.  Nevertheless,  some  of 
these  microphotographs  allowed  us  to  estimate  the  mentioned 
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approximate  size  of  the  metallic  particles  forming  the  colloid. 

For  the  final  solutions,  we  used  the  mixtures  listed  in  the 
following  Table  IA  that  lists  concentrations,  type  of  solvents 
used,  and  approximate  size  of  the  particles  obtained  from  the 
vibrating  orifice  generator  after  evaporation  of  the  dilutants. 

The  size  of  the  particles  listed  in  the  table  has  been 
obtained  estimating,  whenever  possible,  the  position  of  the  first 
minimum  in  the  measured  forward  light  scattering  and  applying  the 
formula  derived  from  the  diffraction  theory  of  light,  x  = 


Mixture  Type 


Solvents:  50*  Water 

50*  Ethanol 

Solution 

Sodium  Chloride* 


Nigrosine  Dye  (75*) 
Ammonium  Sulfate  (25*) 
Solvents:  75*  Methanol 
25*  Water 


Chlorauric  Acid/Sodium 
Citrate 

Standard  Solution 
Solvents:  50*  Ethanol 
50*  Water 


Chlorauric  Acid 
Citric  Acid 
Solvents:  50*  Ethanol 
50*  Water 


Silver  Nitrate 
Sodium  Borohydrade 
Solvents:  50*  Ethanol 
50*  Water 


TABLE  IA 


Reaction 


Concentration 


Fig  Number 
(Appendix  C 


N/A 

6.7xl0~4 

12.6 

la 

6.7xl0~4 

11 

lb 

5.3xl0~5 

7.8 

lc 

6.3xl0~5 

5 

Id 

N/A 

Co=8.5xl0  4 

8.7 

3a 

Co/2 

8.2 

3b 

Co/4 

8.  1 

3c 

Co/8 

7.9 

3d 

Co/16 

6.2 

3e 

Co/64 

6.4 

3f 

Co/128 

3.9 

3g 

Ref  3s 

Cos5.5xl0  4 

Co/4 

9 

4a 

Co/8 

6.7 

4b 

Co/16 

6.7 

4c 

„  -4 

Ref  3a 

Co=4 .3x10 

Co/2 

5a 

Co/4 

5b 

Ref  3a 

Co=8xl0  4 

Co/2 

13 

6a 

Co/4 

11 

6b 

Co/8 

- 

6c 

Co/8 

9 

6d 

Co/32 

6 

6e 

Co/64 

5.4 

6f 

h 


♦Measurements  performed  during  a  previous  research  and  reported  here 
for  the  analysis  in  Legendre  polynomial  coefficients. 
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FISHEYE  POLAR  NEPHELOMETER 


During  the  last  few  years,  various  designs  for  instruments 


capable  of  simultaneously  measuring  the  light  scattering  of 


aerosols  at  different  angles  have  been  proposed.  Reported  in 


this  note  are  some  data  relative  to  a  preliminary  study  performed 


on  a  specific  kind  of  such  an  instrument  that  uses  a  fisheye  lens 


as  a  basic  optical  element. 


Measurements  made  comparing  a  prototype  of  this  instrument 


with  another  instrument  already  tested  will  be  discussed  on  the 


basis  of  the  characteristics  of  the  lens  given  by  the 


manufacturer.  The  performed  observations  have  shown  that 


reliable  measurements  of  light  scattering  can  be  made  in  the 


interval  of  angles  between  5  and  175°,  and  that  the 


characteristics  of  the  lens,  as  given  by  the  manufacturer,  appear 


not  to  be  accurate  at  angles  close  to  45  and  135". 


The  design  of  the  tested  nephelometer  (developed  by  Dr.  G. 


W.  Grams  of  the  School  of  Geophysical  Sciences  at  Georgia  Tech) 


is  depicted  in  Figure  la.  Figure  lb  shows  schematically  the 


operating  principle  of  the  fisheye  lens  prototype  of  the 


nephelometer  herein  discussed,  and  Figure  2  shows  the  geometry  of 


the  system.  This  nephelometer  can  properly  operate  if  the 


particles  are  uniformly  mixed  in  the  sample  volume,  and  if  they 


are  at  a  concentration  low  enough  to  disregard  the  effect  of 


multiple  light  scattering.  Good  uniform  mixing  can  be  obtained 


by  reducing  the  sample  volume  to  limits  compatible  with  the 
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characteristics  of  the  fisheye  lens  employed. 

The  following  discussion  pertains  to  the  specific  lens  used 
during  the  measurements  (Nikkor  Orthographic  Projection  10  mm 
f5.6);  nevertheless,  similar  discussions  can  be  easily  extended 
on  other  cases.  According  to  manufacturer  specifications  on  the 
lens,  the  lens  accommodates  the  180°  of  the  field  of  view 
following  a  formula  that,  for  the  angle  <t>  =  0,  becomes  (see  Fig. 


x  =  L  sin  £ 


(x  in  millimeters; 
L  constant  of  the 
lens  -  10  mm) 


( 1  •  B ) 


where  x  is  the  distance  from  the  intersection  of  the  optical  axis 
of  the  lens  with  the  image  focal  plane  (center  of  the  image)  and 
£  is  the  angle  with  the  optical  axis.  It  is  also  stated  that  the 
radiances  read  on  the  film  are  proportional  to  the  cosine  of 
incidence.  Therefore,  by  changing  the  angle  from  £  to  the  light 
scattering  angle  9,  and  by  taking  into  account  the  geometry  of 
the  laser  beam  with  8,  it  follows  that  the  lens  should  have  the 
following  properties: 

1.  The  radiances  measured  are  directly  proportional  to  the 
particle  phase  functions.  In  fact,  the  radiance  R, 
received  in  front  of  the  lens  in  direction  8,  is: 


R(8,  4>  =  0)  =  o(8)/sin  8 


( 2  .  B ) 


After  the  lens, 


SseS 


ff 


«  y"  . 


R(x,  y 


0)  <r  a  ( 0 )  cos  C/sin  0  =  a  (0) 


(  3  .  B  ) 


The  light  scattering  cross  sections  of  the  particles  are 


directly  proportional  to  line  integrals  of  the  measured 


radiances.  This  is  because  dx  =  d(cos0),  then: 


®  (  0)  d  (  COS0)  a 


•-L 

R(x; 


(4.B) 


The  half  ellipsis  that  results  on  the  image  plane  by 


observing  a  laser  beam  with  the  optical  axis  orthogonal 


to  the  lens  -  but  non  crossing  the  laser  beam,  as  in 


Figure  3 


once  projected  on  the  major  axis  of  the 


ellipsis  will  again  follow  the  same  functional 


dependence.  The  projection  formula  can  be  easily 


visualized  as  projecting  maximum  circles  of  a  sphere  on 


a  plane;  it  follows  that 


x  =  L  cos 0 


y  =  -L  sin0  sin$ 


( 5 .  B ) 


As  further  consequence,  we  can  calculate  the  angular  resolution 


obtainable  from  the  lens  on  the  image  plane: 


dl/d0  =  (d  x2  +  d  y2 )  /d0 


( 6 .  B  ) 


=  L(sin20  +  sin2$cos20) 


For  small  values  of  ♦  : 


dl/d0  »  L  (  s  in2  0+4>2  cos2  0) 


(7.B) 


That  is,  the  resolution  has  a  maximum  for  0  =  90°,  and  is  small 
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the  optical  axis  of  the  lens  and  the  laser  beam  are  perpendicular 
( <t>  =  0°),  dl/de  =  0  for  0  =  0  and  180°.  Therefore,  to  obtain  a 
good  resolution  at  low  and  high  light  scattering  angles,  it  is 
necessary  to  observe  the  image  of  the  laser  beam  with  some 
elliptical  deformation  on  the  x-y  plane. 

Measurements  made  with  a  prototype  of  the  fisheye 
nephelometer  have  been  carried  out  using  a  photographic  camera 
and  a  microdensitometer  for  analyzing  the  films.  A  TSI 
monodisperse  aerosol  generator  has  been  used  for  intercalibrating 
the  prototype  with  the  polar  nephelometer.  The  photographic 
density  (in  the  linear  region  of  the  response)  is  related  to  the 
received  radiances  by  the  relation: 


R ( x , y )  =  a10rD(x'y) 


( 8 .  B  ) 


where  a  and  y  are  constants  that  depend  respectively  on  the 
exposure  and  on  the  type  of  film,  and  D(x,y)  is  the  film  density. 

We  can  observe  now  that  by  repeatedly  scanning  the  image 
(along  parallel  lines  very  close  to  each  other  in  any  convenient 
direction),  it  is  possible  to  correct  the  effect  introduced  by 
the  angle  4>.  In  particular,  if  each  scan  is  the  plot  of  the 
photographic  density  along  segments  (xi,  y^  -  X2,y^,  i  =  l,n). 
The  envelope  of  the  maxima  of  the  measured  densities  is  directly 
related  to  the  light  scattering  through  an  equation  of  the  kind 
of  equation  3  because  the  maxima  are  the  only  positions  where  the 
exposed  portion  on  the  film  completely  intersects  the  probing 
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beam  of  the  densitometer.  Conversely,  a  generic  curve  cannot  be 
correlated  with  the  light  scattering  properties  of  the  particles 
unless  the  value  of  angle  is  exactly  known  ard,  therefore,  the 
response  as  a  function  of  position. 

The  results  obtained  from  numerous  dens i t omet r i c  scans  in 
one  of  the  measurements  is  reproduced  in  Figure  4.  These 
measurements  have  been  obtained  with  monodisperse  aerosol  of 
spherical  particles  of  about  0.8  jjd  radius  (size  parameter  about 
7.9)  with  75%  Ammonium  Sulfate  and  25%  Methylene  Blue.  The 
densitometer  was  set  for  a  resolution  of  about  0.06  pm. 

Figure  5  shows  the  curves  representing  the  envelope  of  the 
maxima  of  the  various  curves,  while  in  Figure  6  the  densitometric 
data  (triangles)  have  been  digitized  and  reduced  following 
Equations  (5)  and  (8).  The  photographic  constants  a  and  r  have 
been  set  to  fit  approximately  the  results  of  the  polar 
nephelometer  (squares).  The  origin  of  the  apparent  differences 
in  the  position  of  the  peaks  has  not  yet  been  fully  understood. 
They  could  be  due  to  variation  in  the  performance  of  the 
viorating  orifice  generator  or  to  non  uniform  distribution  of 
particles  density  within  the  sampling  volume  (contained  in  a  box 
about  4  m  long  and  50  cm  x  30  cm  wide). 

To  assess  the  angular  performance  of  the  lens,  we  have 
performed  measurements  with  OOP  particles  of  size  parameter  about 
12  for  the  two  orthogonal  planes  of  polarization.  Figure  7  shows 
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picture  taken  with  the  fisheye.  Taking  into  account  the  fact 
that  the  polar  nephelometer  cannot  resolve  peaks  and  valleys  in 
the  phase  functions  with  a  precision  better  than  5°,  we  can 
observe  that  there  is  a  slight  systematic  tendency  of  the  fisheye 
coordinates  to  stay  below  the  45°  correlation  line,  especially 
for  directions  around  45  and  135°. 

We  can  conclude  that  the  photographic  measurements 
constitute  a  preliminary  verification,  made  with  a  rather 
primitive  experimental  set  up,  that  a  fisheye  polar  nephelometer 
can  be  successfully  realized  on  the  basis  of  the  observations  and 
discussions  herein  performed.  In  the  present  configuration  the 
system  resulted  to  easily  provide  reliable  measurements  of  the 
aerosol  scattered  light  at  scattering  angles  of  5  and  175°,  with 
an  improvement  of  5°  over  the  laboratory  instrument.  This 
feature  definitively  proves  that  the  development  of  an  actual 
instrument  for  light  scattering  measurements  based  on  this 
principle  is  worthwhile  to  be  pursued.  The  optical  properties  of 
the  lens,  as  given  by  the  manufacturer  can  be  considered 
acceptable  in  the  first  approximation.  A  calibration  should  be 
performed  on  an  actual  system. 
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Interconparison  of  the  results  obtained  with  the  two 
nepheloneters .  (a)  results  from  the  nepheloaeter  in 
Figure  1(a),  (a)  results  fron  the  fisheye  nepheloaeter 
(Figure  1(b)). 
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APPENDIX  C 


MEASURED  PHASE  FUNCTIONS 

The  angular  light  scattering  measurements  herein  reported 
have  been  corrected  for  the  effect  of  the  scattering  angle  on  the 
field  of  view  of  the  detector,  and  for  the  molecular  contribution 
to  the  light  scattering. 

We  have  normalized  the  results  of  angular  light  scattering 
dividing  by  the  best  possible  estimate  of  the  total  light  energy 
scattered  at  all  angles.  In  fact  the  number  density  of  the 
particles  flowing  in  the  sample  volume  depends  on  experimental 
parameters  (such  as  the  performance  of  the  particle  generator  and 
the  amount  of  air  used  for  dispersing  the  particles)  that  can  not 
be  measured  with  sufficient  precision  to  permit  dependable 
determination  of  the  particle  cross  sectional  area.  Therefore, 
we  have  simply  divided  our  angular  data  set  by  the  quantity: 

N 

a=gQ=  I  wsi  V**  (N=33 ;  p=l . 2)  (l.C 
i 

S  p 

where  w  i  are  suitable  weights,  the  measured  values  of  light 
scattering,  and  p  (=1,2)  the  index  for  parallel  and  perpendicular 
polarization  of  the  incident  laser  light. 

In  order  to  compute  the  Legendre  polynomial  coefficients  we 
have  to  take  into  account  that  the  Legendre  polynomial  are  not  an 
orthogonal  set  of  polynomials  in  the  discrete  space  of  our 
experimental  data  (33  points  taken  every  5°  in  the  interval 
between  10  and  170°).  That  is  the  quantities: 


"ij=  I  (i  =  l»n;  j  =  l ,  n;  N=33;  pk  =  cos  0k ) 

k1 


( 2 .  C 


are  different  from  the  analytical  values  of  the  corresponding 
integrals: 


IT.  .  =  P.P  .( 

1J  J,  1  J 


( 3 .  C 


To  obtain  the  best  estimate  of  the  Legendre  polynomial 

coefficients  from  our  data,  we  have  defined  the  matrix  IT,  with 

dimensions  nx33  and  with  elements,  n. .,  so  that  the  equation 

^  J 


g=iinTnii_1uniiTv 


(4 .  C 


provides  the  desired  results  in  the  least  square  sense. 

In  this  study  the  optimization  of  the  series  has  been 
performed  with  a  trial  and  error  method  by  varying  the  dimension 
n  of  the  array  of  the  coefficients,  g,  and  of  the  matrix  JI.  As 
optimum  value  of  n  has  been  taken  the  one  whose  associated 
polynomial  was  best  fitting  the  measured  values  without 
exhibiting  physically  irrisonable  oscillations  in  between  the 
experimental  points.  Naturally,  varying  the  order  n  of  the 
polynomial,  different  values  for  the  low  order  coefficients  are 
found.  In  particular  the  value  of  the  coefficient  gQ,  varies 
from  the  expected  unitary  value.  This  basically  indicates  the 
existing  difference  between  the  normalization  constants  of  the 
experimental  data  and  of  the  optimized  polynomial  series. 
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Linear  inversions  performed  with  eq  4.c  are  well  known  to 
produce  negative  and  physically  meaningless  solutions  when 
applied  to  experimental  data  affected  by  experimental  errors  or 
by  insufficient  angular  resolution.  Therefore,  unconstrained 
analysis  of  this  kind,  when  giving  negative  solutions,  will 
provide  a  necessary  (but  not  sufficient)  indication  of 
applications  on  ill  conditioned  cases. 

The  figures  that  follow  report  the  angular  light  scattering 
measurements  for  the  various  aerosols  studied.  The  triangles  and 
the  circles  in  the  figures  are  respectively  the  experimental 
points  for  polarization  of  the  incident  light  parallel  and 
perpendicular  to  the  plane  of  scattering  (the  crosses  have  no 
scientific  content  since  are  numbers  proportional  to  the  recovery 
time  of  the  photomultiplier).  The  continuous  lines  are  the  fits 
of  the  experimental  data  computed  with  the  Legendre  polynomials. 
The  set  of  Legendre  polynomial  coefficients  for  the  two  curves 
are  reported  underneath  each  figure  (first  column  the  fit  for 
triangles,  second  for  circles).  Some  of  the  times,  a  third  column 
is  present  with  the  averages  between  the  first  two  set  of 
coefficients  (unpolarized  light). 

In  general,  by  observing  the  figures  it  is  possible  to 
notice  that  the  computed  polynomials  fit  the  experimental  data 
better  in  the  case  where  the  experimental  phase  functions  have 
oscillations  of  small  amplitude.  This  occurs  in  general  in  the 
case  of  small  particles  or  in  the  case  of  particles  with  complex 


internal  structure  or  shape. 

In  the  case  of  small  particles  this  is  not  surprising,  and 
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the  reasons  have  been  discussed  in  page  7.  In  the  case  of  large 
particles  with  a  complex  shape  or  internal  structure,  the  results 
can  be  interpreted  as  due  to  the  combined  effect  of  the  low 
angular  resolution  of  the  experiments  and  of  the  random 
orientation  of  the  particles.  Therefore,  in  our  opinion,  the 
reported  polynomial  coefficients  in  the  case  of  these  particles 
cannot  be  considered  any  longer  an  experimental  estimate  of  a  set 
of  true  coefficients  that  can  be  uniquely  associated  with  the 
particles  under  study.  Sometimes  the  polynomials  fitting  the 
measurements  for  the  two  polarizations  of  the  incident  light  do 
not  merge,  as  they  should,  at  0  and  180°.  At  this  stage  of  the 
research  we  have  performed  separate  fits  for  the  two  sets  to 
better  the  actual  applicability  and  versatility  of  the  method. 
In  fact,  we  did  not  include  in  our  equations  the  fact  that  the 
two  set  of  data  are  actually  describing  the  same  phase  function, 
a  feature  that  can  be  easily  taken  care  by  including  in  the 
analytical  formulae  the  azimuthal  dependence  of  the  scattered 


light . 

Nevertheless,  this  method  has  been  chosen  to  show 
independently  for  the  two  polarizations  of  the  incident  light  the 
few  cases  where  the  linear  inversion  method  yields  negative 
solutions  (fig.  3a,  c,  e,  f;  5a,  b;  6a).  As  it  is  well  known 
from  the  Lorenz  Mie  theory,  this  fact  can  be  expected  to  occur 
more  frequently  in  the  case  of  large  particles  or  of  relatively 
high  light  absorbing  particles  because  they  are  those  that 
possibly  require  more  terms  in  the  Legendre  polynomial  series. 
However  it  can  be  noticed  how  negative  coefficients  never  appear 
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FIGURE  1 
Sodium  Chloride 
(50%  Ethanol-50%  Water) 
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FIGURE  3 

Nigrosine  Dye  (75%) -Ammonium  Sulfate  (25%) 
(75%  Methanol-25%  Water) 
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FIGURE  4 

Chlorauric  Acid-Sodium  Citrate 
(Prepared  as  in  Ref  3A  and  diluted 
in  50%  Ethanol-50*  Water) 
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FIGURE  5 

Chlorauric  Acid-Citric  Acid 
(Prepared  as  in  Ref  3A  and  diluted 
in  50%  Ethanol-50%  Water) 
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FIGURE  6 

Silver  Nitrate-Sodium  Borohydrade 
(Prepared  as  in  Ref  2A  and  diluted 
in  50*  Ethanol-50*  Water) 
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